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TECHNICAL NOTE 2875 


BEHAVIOR IN PURE BENDING OF A LONG MONOCOQUE 
BEAM OF CIRCULAR-ARC CROSS SECTION 


By Robert W. Fralich, J. Mayers, and Eric Reissner 
SUMMARY 


Án analysis is made of the behavior under a loading of pure bending 
moment of a thin, infinitely long, pure-monocoque beam having a constant, 
doubly symmetric, circular-arc cross section. Bending moments, deflec- 
tions, and stresses are obtained. The analysis shows a nonlinear 
behavior in bending which leads ultimately to a maximum moment and insta- 
bility. 


INTRODUCTION 


The pure-monocoque beam, with no internal bulkheads, ribs, spars, 
or stiffeners, represents a limiting structure which designers can 
approach in an attempt to obtain thin hollow wings with low fabrication 
and assembly costs. Examination of the structural behavior of the pure- 
monocoque beam is, therefore, of interest in order to determine its 
favorable or unfavorable characteristics. In particular, an understanding 
of the nature and growth of the cross-sectional distortion that arises 
from the tendency for the beam to flatten under longitudinal bending 
loads is important. With sufficient flattening, the beam bending stiff- 
ness can be reduced to a point at which the beam no longer will sustain 
an increase in bending moment, and instability results. The related 
problem in connection with circular cylinders is well-known and has been 
discussed by Brazier (ref. 1). 


In the present paper a thin, infinitely long, pure-monocoque beam 
having a constant, doubly symmetric, circular-arc cross section is ana- 
lyzed under a loading of pure bending moment (fig. 1). This loading 
produces identical deformation of all cross sections and, therefore, 
permits a relatively simple structural analysis to be made. Elastic 
behavior is assumed and local buckling is not considered. 


ln an actual pure-monocoque circular-arc wing, these idealized con- 
ditions of structure and loading would not be realized - the cross sec- 
tion and bending moment would vary spanwise and the bending moments 
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would be produced by lateral forces applied directly to the wing surface. 
At least the order of magnitude of the distortion of any cross section 
in the actual wing (away from the root or tip bulkheads) could probably 
be determined, however, by assuming the cross section to be part of an 
infinitely long uniform beam subjected to a uniform bending moment equal 
to the actual local bending moment. 


An analysis of the idealized structure is given in the appendix. 
The body of the paper contains the results of the analysis and describes 
the distortion characteristics and stresses of the pure-monocoque beam 
under uniform bending moment. 


Part of the work presented herein was carried out in June 1951 
while the last-named author was temporarily at the Langley Laboratory, 
and it was continued by correspondence. 


SYMBOLS 

A, B, C constants 
b semiwidth of cover plates, in. 

Et3 
D flexural stiffness of cover plates, A, ov in-1b 

12(1 - pê) 

E Young's modulus for beam material, 1b/in.* 
h rise of circular arc measured from undeformed middle sur- 


face of beam to middle surface of covers, in. 


Lo moment of inertia of undeformed beam cross section, in." 

k curvature of beam axis, in zl 

My; M, plate bending moments on cross sections perpendicular to x- 
and y-axes, respectively, in-1b/in. 

m total beam bending moment, in-lb 

Ny» Ny plate middle-surface forces in x- and y-directions, lb/in. 


Qy plate shearing force in xz-plane, 1b/in. 


= maim eee eee ze 
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Subscripts: 
b 


t 


radius of curvature of circular-arc cross section, in. 
thickness of cover plate, in. 


radial deflection of point in middle surface of cover of 
bent beam, in. l 


total radial deflection of point in middle surface of cover 
plate, in. 


coordinates 
parameter, y = Ab 


plate axial strain 


strain at surface passing through beam edges 


parameter defined by KE BUE. in. ^l 


Poisson's ratio for beam material, wp = 0.316 for all 
computations 


local buckling stress of compression cover 


stress functions 


linear differential operator defined by 
T n h 
ie e Rie 
ax* — Aë oy! 


bottom or tension cover 


top or compression cover 


oo ae e. ee een osa 
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RESULTS 


Relationship of Bending Moment and Beam Curvature 


For calculations of the over-all stiffness of a beam, a knowledge 
of the relationship between bending moment and longitudinal beam curva- 
ture is important. This relationship for the pure-monocoque beam is 
shown nondimensionally in figure 2. 


Unlike solid beams or shell beams with internal stiffening, which 
have linear moment-curvature relatiqnships, the pure-monocoque beam 
exhibits a nonlinear relationship. This nonlinear behavior is a conse- 
quence of the continual reduction in cross-sectional moment of inertia 
as bending progresses. This reduction in moment of inertia or resist- 
ance to bending leads ultimately to & maximum moment. The slight spread 
in the vicinity of the maximum among the curves for different values of 
t/h is a result of taking into account the shift in the neutral surface 
of the beam due to the unsymmetric behavior of the top and bottom covers. 
If this small effect is neglected, the maximum value of the bending 


Etlo 
rh ? 


moment is m = 0.285 corresponding to a nominal midchord critical 


stress of 0.285 = based on the geometry of the undeformed cross sec- 


tion, or a midchord critical stress of 0.36 9» based on the geometry 
of the flattened configuration. This maximum value of bending moment 
which occurs at a beam curvature of k = 0.58 + may be considered to 
be the upper limit of the strength of the beam. 


Distortion of Cross Section 


A measure of the cross-sectional distortion is given by the magni- 
tude of the beam-thickness change at midchord. The ratio of midchord 
beam thickness of the deformed cross section h4(0) + h,(0) to midchord 


beam thickness of the undeformed cross section @h is plotted in fig- 
ure 3 a8 & function of the beam-curvature parameter. From the figure, 
the flattening is seen to be appreciable over most of the loading range, 
and at the value of the beam-curvature parameter corresponding to the 
maximum bending moment, flattening has reduced the original beam thick- 
ness by almost 30 percent. The effect of different t/h ratios does 
not appear in this plot since the thickness change is independent of 
neutral-surface shift. 
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The cover deflections at every chordwise point may also be of 
interest. These deflections can be computed from equations (48) and (49) 
or equation (65) in the appendix. A typical chordwise plot of the cover 
deflections, in which the small effect of the neutral-surface shift is 
neglected, is shown in figure 4 for a beam loaded very nearly to its 
maximum moment. 


Stresses 


The magnitudes of the axial stresses in the beam are best indicated 
by the magnitudes of the axial stresses at the midchord where they have 
their greatest values (at least up to maximum moment). For each cover, 
the ratio of midchord axial stress to the linear-theory stress at the 
same beam curvature is plotted in figure > as a function of the nondimen- 
sional beam-curvature parameter for values of t/h equal to 0 and 0.4. 
When compared at equal beam curvatures, the stresses of the present the- 
ory are always of less magnitude as a consequence of the flattening. 


The ratio of midchord axial stress for each cover (with the small 
effect of neutral-surface shift neglected) to the linear-theory stress 
at the same bending moment is plotted in figure 6 as a function of the 
nondimensional bending-moment parameter. The result indicates that, at 
the same value of bending moment, the cover stress in the present theory 
is greater than the linear-theory stress. This condition is a result of 
flattening decreasing the moment of inertia more rapidly than it decreases 
the beam thickness. 


DISCUSSION 


The results indicate a nonlinear relationship between bending 
moment and beam curvature that arises from the low chordwise bending 
stiffness of the covers and ultimately leads to an instability. The 
results also indicate that, in general, an appreciable profile flexi- 
bility exists which may be detrimental from an aerodynamic point of view. 
In order for completely hollow wings to resist large distortions, addi- 
tional chordwise bending stiffness must be realized in the covers (e.g., 
by integrally stiffened or sandwich construction). 


In the analysis the effect of chordwise middle-surface stresses on 
the deflections is neglected. Order-of-magnitude considerations indi- 
cate that such an assumption is justified for cross sections where 


(b/r)? and (t/n)? are smali compared with unity. These conditions 
are satisfied, in general, for shallow sections such as would be con- 
sidered for thin wings. 


6 | NACA TN 2815 


The effect of local buckling of the compression cover has been 
omitted from the analysis. Its importance with respect to the present 
problem can, however, be estimated in the following manner: 


The theoretical buckling stress in uniform compression of & long 
curved plate is given by the relation 


_ Et 
SX = 0.6 s 


which applies to both simply supported and clamped long plates if the 


2 
lib V 2 
curvature parameter = V1 - is eater than about 60. In this 
p ¡ u gr 


range more than one buckle appears in the circumferential direction and 
the buckling stress becomes practically independent of edge conditions. 


It is well-known, however, that the "classical" stress of 0.6 = is 


never attained in actual curved plates and requires modification by an 
appropriate "knockdown" factor. This factor is a function of the ratio 
of radius to plate thickness and has been discussed in many papers, for 
example, reference 2. In this reference the buckling stress is given by 
the empirical relation 


me. = (0.588 - 0.000433 t) = 


for plates having radius-thickness ratios between 500 and 1,000. The 
knockdown factor in the region r/t > 1,000 can be taken as the same or 


less than that at = = 1,000. 

Local buckling of the compression cover in bending is estimated to 
occur if the stress measured at the midchord is between 20 and ho per- 
cent higher than the stress required for local buckling in uniform com- 
pression. The estimates are based on curved-plate linear-theory con- 
siderations and experimental results for cylinders (ref. 3). If an 
average value of 1.3 is assumed for the ratio of the critical stresses, 
the stress for local buckling in bending becomes 


rl Et 
Ox p = (0.765 - 0.000563 z) P 


= ——. — t e e 
o — —ÀMÀ—— a a á———— I 
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Since the midchord stress corresponding to flattening instability is 


approximately 0.36 = (gee fig. 5), local buckling of the compression 


cover would be expected to occur after flattening instability for 

r/t <720 and to precede flattening instability for r/t > 720. By 
taking flattening into account, however, this value of r/t is decreased, 
since the chordwise curvature of the covers, and hence, the critical 
compressive stress, is decreased. 


CONCLUDING REMARKS 


An analysis is made of the behavior under a loading of pure bending 
moment of a thin, infinitely long, pure-monocoque beam having a constant, 
doubly symmetric, circular-arc cross section. Bending moments, deflec- 
tions, and stresses are obtained. The analysis shows a nonlinear rela- 
tionship between bending moment and beam curvature that leads ultimately 
to a maximum moment and instability. If no interaction is assumed between 
local buckling of the compression cover and flattening instability, local 
buckling should precede flattening instability of the beam when the 
radius-thickness ratio of the covers is greater than about 720 and flat- 
tening instability should precede local buckling when this ratio is 
smaller than about 720. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., October 16, 1952. 
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APPENDIX 


ANALYSIS OF IDEALIZED STRUCTURE 
UNDER LOADING OF PURE BENDING MOMENT 


Statement of Problem and Assumptions 


The behavior of the idealized structure under a loading of uniform 
bending moment can be determined from consideration of the equilibrium 
conditions for a small element in a chordwise strip in each of the covers 
of the bent beam (see fig. 7). A reduction of the problem to the analysis 
of a chordwise strip is justified on the basis of symmetry considera- 
tions, since all chordwise sections behave identically under the given 
loading conditions. The small elements and the resultant forces and 
moments acting upon them, together with the sign convention used, are 
shown in figure 7. The circumferential middle-surface forces Be, 


and Nyy are not included in the analysis since order-of-magnitude cal- 


culations indicate that, for shallow cross sections, the effect of the 
circumferential forces on the deflections are not of first-order impor- 
tance. 


The relative simplicity of the following analysis and the resulting 
explicit formulas, in contrast with Brazier's work on the circular cylin- 
drical shell, are due to the fact that the shell sections considered are 
shallow. The present analysis may be readily extended to other than 
circular-arc sections by retaining the shallow-section assumption. 


Derivation of Basic Equations 


Since all chordwise sections behave identically under the given 
loading condition, only force equilibrium in the z-direction and moment 
equilibrium about the x-axis have to be considered. For a small element 
in a chordwise strip of the top or compression cover, these equilibrium 
conditions are 


t 
A Nx, k = 0 (1) 


2Y 
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“yt 


dy ` Kéi adis (2) 


where -k is the beam curvature, which is constant along the length of 
the beam. For a small element in a chordwise strip of the bottom or 
tension cover, the equilibrium conditions are 


y 7 Mek = 0 (3) 
and 
a db. (4) 


Elimination of ba and Ho, from each of the pairs of equations yields 


dEM 
Yt 
ay? - Nx,k = 0 (5) 
and 
d^ My, 
WB ako (6) 


Since all chordwise sections behave identically, the cross-sectional 
distortion is independent of the axial coordinate x. As a result, for 
each of the chordwise strips, the moment-distortion relationships from 
plate theory become 


aM GL 


My, Se (7) 


— ES eto ee eiie ak ri dei — o — — ——————— wf A ———— S — = aa aana a — —À E ge ee i E 
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Hy = EO p (8) 


+ 


where the radial deflections Wi(y) and Wb(y) are positive when 
directed downward in figure 7. Substitution of these relations for My 


and My. into the two equilibrium equations (5) and (6) yields 


a'W.(y) — _ 
-D De Ny k = 0 (9) 
dy 
and 
-D a y) - Ny, k = 0 (10) 
ay 


The axial forces EA and Ba, can be written in terms of the 


axial strains as 


Ny, = Etex, (11) 
and 
Ta = BYE x, (12) 


The strains, measured from the as yet undetermined neutral surface of 
the cross section, are 


€. = zl E 3 - x + kW, (y) (13) 
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and 
E, = a - a - Ex + ki, (y) (14). 
where € is the strain at the surface passing through the beam edges 


Xo 


and dé - a is the distance from this surface to the middle surface 
b 


of the undeformed covers. The latter quantity gives the rise at any 
point y in good approximation for shallow sections. With the fore- 
going relations, the equilibrium equations (9) and (10) become 


d sU) y^ 
-D "ach + Etk de - a + Exo - KW. (y) = 0 (15) 
and 
-D a Ww) - Etk gd - 5) - Ex, + kWy(y)} = 0 (16) 
ay v^ 
Rearranging the terms gives 
am (y) y2 
p E tn, (y) = PAI - 55) + Etke (17) 
SEN 2 Xo 
y b 
and 
a'w y) 2 
D B + Etk m (y) = -Etk^h l- 5) + Etk€x (18) 
dy pe O 


The two foregoing equations are the differential equations of equi- 
librium for the chordwise strips in each of the covers. Each equation 
may be noted to be that of a beam of unit width on an elastic found&- 
tion, where the left-hand side of each equation represents the restoring 


= .——-— — MÀ mee — —á— ` e ~ i a A A A ATT = AA me re M 
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forces due to the cover bending stiffness and the elastic foundation 
and the right-hand side represents the forces causing deflection. 


The deflections can be obtained from the equilibrium equations for 
a given set of boundary conditions, which in this problem are taken to 
be 


GÉIE) EC O 
AW¿(+b) aw, (tb) 
E GN 
and 


a&W,(3b) —a*W,(sb) 
—————— + M 


ay? dy? 


These conditions indicate that no edge deflections exist relative to the 
bent beam, that the edge angles are maintained between the two covers, 

and that no resultant edge moments are permitted. The last condition is 
obtained from My, (+b) + My, (+b) = O by virtue of the moment-distortion 


-2uk + 


relationships. 


Before the deflections can be completely defined, the strain at the 
surface passing through the beam edges is required. Because of symmetry, 
no resultant axial force exists in the beam at any cross section. This 
condition is given by the relation 


b D 
I, Ny. dy dp Ny, dy = 0 (19) 


and leads to the relation for the strain at the surface passing through 
the beam edges 


b 
Ex, = ši |. + Wily)| dy (20) 


when the force-strain relations are substituted into equation (19). 
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Alternate Derivation of Basic Equations 


The basic differential equations can be derived by considering the 
total deflection of small elements in chordwise strips of each of the 
covers measured from the initially nonstressed position of the elements. 
Since the deflections measured in this manner will be large compared 
with the cover thicknesses, the equilibrium conditions for the small 
elements must be expressed by the equations of large-deflection theory 
of shallow cylindrical shells (ref. 4). 


For the top cover, these equations are 


mary) 1% , 3 n. Guy) [8214 ( y) 


l va (x,y) = —— 
ST Ss a e 372 


+2 (21) 


and 


34g. (x, y) Ee (x,y) " py (x, y) Pw. (x,y) adj 


d 
ay? sa ee ay? 
"T ^g. (x,y) 9%, (x, y) 


D huel e E 


dx dy dx dy (22) 
and for the bottom cover, the equations are 
2 
M NEM Im (xy) |" Ow (xy) d (x,y) E (23) 
pU pm Ox Oy Ixe dy? zi 
and 
li _ pp (x,y) 32m, (x,y) dpp (x, y) [dm (x,y) 1 
Dvw(xy) =t 52 52 EuS ae "€ - =| - 
2 
T d Py (x,y) dw (x, y) (24) 


Ox dy Ox dy 


14 
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where wy (x,y) and wy (x, y) are the total deflections and plx y) 


and Pp (x, y) are the respective stress functions of the covers. 


The 


relationships between the stress functions and the middle-surface forces 


in the covers are 


pj (xy) B "Yt Ip (x, y) _ To IE (x, y) B Tore 
óc ° aye t E 
(25) 
Ip, (x, y) _ "Yb Ip, (x,y) _ KI op, (x, y) _ Ten, 
SN ày? t dx dy t 


These equations together with the boundary conditions determine 


W¿(x,y), Wip (2, y); Qux y), and Pp (x,y). 


If use is made of the iden- 


tical behavior of each cross section, the deflections of the elements 


may be assumed to be 


we(x,y) = -5 ka + Wely) (26) 
and 

mo (x,y) =-5 kx^ + Walz? (27) 
where d ka is the parabolic deflection of the beam edges due to the 


uniform bending moment and W,(y) and W,(y) 
tions of the covers. Since the circumferential forces 


are not included in the analysis, the stress functions 


g(x,y) = (y) 


Pp (x,y) SV (y) 


are the chordwise deforma- 


nyc and Ny, 


may be taken to be 


(28) 


(29) 
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Substitution of the relations for the deflections and the stress 
functions into the large-deflection equations gives 


à 2y 
ay ay? 
2 
D a Wely) wak ady) (31) 
dy“ dy2 
and 
4 2 
ob LI o 4 bias 
dy ay? 
am (y) ` ay ol) "a" 
ay" ay? 


After the first equation in each case is integrated twice, the result is 


2 
ao, (y) 2 
E Wily) + vo} + Ciy + Co (34) 
E dy? 2r 
and 
2 
aso, (y) 
1 ONY ye 
E = Mwy) - | + cy +c (35) 
E aye E 3 39 + Cy 
From symmetry considerations, Cy = C3 = O so that substitution of these 
2 
ao, (y) a“, (y) 
expressions for = PE into the remaining two equations 


gives 


s ol. er ee ie: a 
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d 
D id + EK Vull = == - EtkCo (36) 
dy 
and 
d 2 
5 SCH + Etk^WL(y) = SA ZS - EEN (37) 


The boundary conditions on the total deflections wy,(x,y) 
and W(X y) are 


Loe 
w,(x,tb) = wp(x,tb) = -5 kx 


Ow, (x,+b) _ ow, (x,+b) 


oy Oy 


and 


My,(x,tb) + My (xb) = 0 


so that the corresponding boundary conditions on W,(y) and W,(y) 
become 


W¿(+b) = M(+b) = 0 


Oe Läb) aw (+b) 
"dx = dy 


2 Boot: 
am, (4b) . a^W, (+b) "^ 


dy? dy? 


-2uk + 
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The conditions necessary for determining Co and Ch are 


Nig, (4b) = Ny, (4b) (38) 


Ip b 
N ay + | Ny, dy = O (39 
-b Xt _b Xp 


Which account for continuity of axial forces at the edges and zero result- 


ant axial thrust. If, in these equations, c and E are written 


in terms of deflections by use of the relationships between axial forces 
and stress functions and between stress functions and deflections, the 
following equations are obtained: 


2 
E. + Cp - Cy = 0 (40) 


" 
JR Gn + mo (y) | dy + 2b (Co + Cy) =0 (41) 


The simultaneous solution of these equations for Co and Cp yields 


b 
Co = -kh - EL Le + we) dy (42) 
and 
b 
Cy = kh - el |? + wo) dy (43) 
: where 
dus BE 
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Substitution of these expressions for Cy and Cj into the equi- 
librium equations (36) and (37) results in 


"n a^i, (y) 


dy 


b 
ača - a + Etk Së us? + Wy(y)| dy 


b 
(4) 


and 


aw, (y) 
dee + Etk (y) e 


b 
Sa - d + Etk el. |n * wy) dy 


(45) 


The two foregoing equations are identical with the equilibrium equa- 
tions derived in the preceding section since the integral appearing in 
these equations is recognized as €, , the strain at the surface which 

O 


passes through the beam edges. 


Solution of Basic Equations 


The differential equations for the deflections of the covers 
(eqs. (17) and (18)) are 


D Dem + Etk Wily) = Etk^h e : d + Bike, 


Y 2 


D SW) Etk^W,(y) = -Etk^h k I 5) + Etke, 
dy b S 
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These equations upon division by D become 


D e 
wt + ha Cy) = d = E) + hat e (46) 
dy b 
and 
amy) h y yê i US 
= +a my) = -ahl -j a > (47) 
dy b 
where 
p, = Etk 
D 
The solutions are 
SE, 20 
Wi (y) = A, cosh Xy cos Ay + By sinh Xy sin Ay + d ~ `z ies 
(48) 
and g 
€ 
Wp (y) = A cosh Ay cos Ay + By sinh Ay Bin Ay - i: - S + E 
(49) 


where only two arbitrary constants are required in each solution because 
of symmetry. 


— A A A A re A ee ee eee — EE —  — ay ee s ef a ==> — 


If use is made of the boundary conditions and the relation for determining €x? the con- 


stants become 


h sinh 27 + gin 2) 
t 
V n sinh 2y + sin 2y)sinh y sin y - (sinh 27 - sin 2y)cosh Y cos (50) 
> (sinh 27 + sin 2y)(cosh 27 + cos 27) 
3(1 - y?) 


(sinh 27 + sin 2y)(cosh 27 + cos 27) 


| - + sin sinh 7 sin ~ (sinh 2y - gin cosh = (51) 
2 


Pt bk cosh 7 cos Y A 
h 7 sinh 27 + sin 27 


sinh - sin sinh gin 7 + (sinh 27 + sin 27)cosh ? cos (52) 
2) (sinh 27 + sin 2y)(cosh 27 + cos 27) 


Oc 


Glgc Mb VOVN 


and 
By 1 _ cosh y cos y a 
h X y Binh 2y + sin 27 
TEC 
h ginh - gin 2y)sinh y sin y + (sinh 27 + sin cosh y cos Y (53) 
e (sinh 27 + sin 2y)(cosh 27 + cos 27) 
3(2 - 12) 
where 
d 2) / xr 
Y = Ab =Xf1211 - — 
( A 
and 
b? = 2rh 
The relation for ex, Obtained from equation (20) and used in evaluating the constants A, 
e = 
Ay, By, and By ie e . HEN eee Se This result shows that the shift of the 


Xo py2 sinh 27 + sin oy" 
neutral surface is in the direction of the tension cover plate, 


With the constants defined, the deflection at any point in the cross section can be deter- 
mined. 


Glge ML VOVN 


Te 
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Beam Bending Moment 


The beam bending moment at any cross section is determined from the 
sum of the moments acting on the cross section. If the bending moments 
of the covers about their own middle surfaces and the bending moments of 


the axial forces Nx, and Ny about the neutral surface of the cross 


section are taken into account, the total bending moment is 


€ 


b b b z 
A IS dy + 


b Ex 
J N. d - d - => + W (y) | ay (54) 


After substitution for Ny, and No as given by the force-strain 


relationships (11) and (12) and for Me, and M, as given by the fol- 
lowing moment-distortion relationships obtained from plate theory, 


A M (55) 
=- pl- 12 
and 
OW 
M E, cid (56) 
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equation (54) becomes 


b 2 
SE » f T mr , d d 


b dy* dy? 
b € i 
Zo 
six | V EL a+ 
b € P 
Bt | d Ej - 2 + W (7) dy (57) 


With the bending moment reduced to a function of the deflections and 
the deflections defined by the equations of the preceding section, the 
bending moment can be evaluated to give in terms of nondimensional para- 
meters 


15(cosh 27 + cos 27) 8 dc 2y - cos 27) | _ 


2 
u 
1673(sinh 2y + sin 27) 8(1 - u2)V/ (cosh 27 + cos 27) 


15(cosh 2y cos 27 + 1) 14 TE BE y-sinh 27 sin 27 
ly*(sinh 27 + sin py)? oh (1 = u? ) DI (cosh 2y cos 27 + 1) 
(58) 
where 
E EN 
t/h 
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and Ip, the moment of inertia of the undeformed cross section, in which 


the moment of inertia of the plates about their own middle surfaces is 
neglected is defined as 


b 2 
NEA 
EEN dE 


mr/EI 
The moment parameter T, plotted against the beam-curvature para- 
t/h 


meter v for several values of t/h, is shown in figure 2. 
t/h 


Cross-Section Distortion 


The magnitude of the cross-section deformation may be given in terms 
of the ratio of beam depth under loading to beam depth of the unloaded 
section measured at the center of the cross section. That is, 


hc(0) + hy(0) _ LÀ sinh y sin 7 


2h u v(sinh 2y + gin 27) (59) 


where 

hi) = d : ei “wis 
and 

h (y) = d - 5) + Wy) 


This ratio, plotted against the beam-curvature parameter, is shown in 
figure 3. 


The shape of the deformed cross section may be found from the equa- 
tions for chordwise deflection of the covers. If the small effect of Ex 
O 


4Y 
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is neglected so that W,(y) = -Wy(y), the deformed shape of each cover 
is given by 


h y(sinh 27+ sin 27) 
(60) 
tI 
This ratio is illustrated in figure 4 for a loading of m = 0.285 a 


Axial Stress 


From a preceding section (eqs. (11) to (14)), the nondimensional 
axial-stress ratios can be written 


N € 
XL y2 Xo Wi(y) 
Eh > |O 32 * "e" h (61) 
and 
NY y? €x. = Wy) 
b _ _ EO b 
= z pp ue (62) 


In general, the nondimensional axial-stress ratios at the greatest 
distance from the neutral surface are of primary interest. At y -90, 
therefore, 


li) 
kor 4 sinh y sin 7 u 


t 
pa - ZG a ee y(coth v - cot y) (63) 
Gi 


- £ y(coth 7- cot 7) (64) 


wl: A u2) 


AE i 
_ o A i re e o i i i —— on = 


hi(y) " hy(y) B (sinh y sin y cosh Ay cos Ay -cosh y cos y sinh Ay sin Ay 
INE ORE 
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The preceding equations are illustrated in figure 5 in terms of 


C) 0 (o) o 

RL and 22 plotted against the beam-curvature parameter 
Etkh Etkh E 

kr _ for several valu T 

t/h es of t/h. 

The curves plotted in figure 5 can be presented in another form by 


referring the stress ratios to the elementary-beam-theory stress at 
equivalent moment rather than at equivalent curvature. Thus, 


(3 y (Ex). kr 
y=0 |. Etkh t/h 
tmh/Io mr/EIQ 
t/h 
and 
("x ) (o) =0 kr 
b/y=0 Etkh t/h 
tmh/Io mr/EIQ 
t/h 


In figure 6, these nondimensional stress ratios are plotted against 
the nondimensional moment parameter, with the small effect of changes in 


the t/h ratio being neglected. 
Approximate Formulas 
When the shift in the neutral surface is neglected, simplified expres- 


sions for deflection, bending moment, and axial stress are obtainable 
from the resulting equations by expanding functions of y into finite 
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power series to give 


Wily) SS Wy) 


h h 
yi y} 
L- 99 poe 2 
SE cosh Ay cos Ay + — — ———— sinh Ay sin Ay + 1 - d 
. + —áÓÀ 
15 15 
(65) 
A wa] LA e" 
1+ E 
15 
and 
-N 
Ec: qu 
Etkh  Etkh 
l - Ce Ts pad 
SN cosh Ay cos Ay + = EEN sinh Ay sin Ay (67) 
1 + az "qu S 
15 15 


The foregoing approximate formulas are considered to be very satis- 
factory in the region 0 < 2:3 $ 0.6. The maximum errors to be expected 


t/h 


are of the order of several percent. 


MM a EE ——— a — IMP ya EQ ga P s PL ILL Le. ——— PM > = e 
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Figure 2.- Relationship between bending moment and longitudinal beam 


2 
curvature. (o = 3e tbh = 8 2) 
15 19 32 


NACA TN 2875 31 


8 
h,(O) +hg(O) 
2h 
6 


Figure 3.- Ratio of beam thickness of loaded section to beam thickness 
of unloaded section. 
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Figure 4.- Chordwise cover deflections under a loading of m = 0.285 ——. 
Neutral-surf&ce shift is neglected. 
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Figure 5.- Ratio of midchord axial stress to elementary-beam-theory 
stress Ekh at equivalent curvature. 
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Figure 6.- Ratio of midchord axial stress to elementary-beam-theory 
stress mh/T, at equivalent moment.  Neutral-surface shift is 
neglected. 
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Figure 7.- Forces and moments acting on cover elements. 
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